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ABSTRACT 



Minimum expected-cost solutions to the single-period 
inventory problem (the Newsboy problem) have been well known 
for many years, for risk cases where the distribution of 
demand is known. Also well known are minimax cost and 
Laplace solutions for the uncertainty case where only the 
range of demand is known. This study explores the use of 
order-statistic-based quantile estimators as decision pro- 
cedures while data gathers during the early decision periods. 
Simulation, using a variety of demand distributions and 
unit-cost values, provides recommendations on which period 
to leave the minimax procedure for a suggested data- based 
decision rule which will be asymtotic to the optimal risk 
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INTRODUCTION 



I . 

The classic single period inventory problem is called by 
various well-known names such as "Newsboy problem," "News- 
paperboy problem," and so on. [6]. In this thesis it will 
be called the Newsboy problem . 

Solutions to the Newsboy problem are well known under 
both "risk" and "uncertainty" conditions. [7]. The purpose 
of this thesis is to study the Newsboy problem under initial 
conditions of "Uncertainty," (i.e., when the possible future 
for demand is not yet known) and then, for the subsequent 
periods as data gathers, to suggest, structure, and evaluate 
decision procedures during this transition phase to decision 
under "Risk." 

The inventory model studied in this thesis is a time 
independent (stationary) and linear-cost type. The study 
will consider the problem of minimizing the expected cost. 
Since with the Newsboy problem both the expected profit 
maximization and expected cost minimization cases have 
identical solutions, the study should be applicable to both 
interpretations of the problem. 

Chapter II reviews the Newsboy problem for decisions 
under risk and decisions under uncertainty, and then suggests 
possible alternative procedures for "from uncertainty to 
risk" conditions. Among these. Order Statistics (a non- 
parametric approach) appear promising and are discussed in 



7 



detail in Chapter III. Chapter IV evaluates the alternative 
nonparametr ic transition phase procedures of the Newsboy 
problem and a$ a result proposes a transition phase pro- 
cedure. In Chapter V conclusions and recommendations for 
further work are given. 
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II . SOLUTION OF THE SINGLE -PERIOD INVENTORY (NEWSBOY) 
PROBLEM UNDER RISK AND UNCERTAINTY CONDITIONS 



In this section the Newsboy problem is examined under 
risk and uncertainty, and optimal decision rules for these 
conditions are given. These two cases form the bounds of 
the transition phase and permit us to characterize the 
transition phase decision problem we intend to study. 

Here, when it is said, "under risk," it is meant that 
the probability distribution of demand is known. When the 
distribution of demand is not known, and data in hand is 
not sufficient to estimate the possible futures of demand, 
then the problem is called "under uncertainty." [7]. 



A. SOLUTION OF THE NEWSBOY PROBLEM UNDER RISK 

Let f^(x) be the probability density function of the 
demand for continuous populations, and p^(x) be the proba- 
bility mass function for discrete populations. If C sQ is 
the stock-out or shortage cost per unit of unsatisfied 
demand, and C the carrying or surplus cost per unit of 
excess supply, then the cost for a particular period of 
Newsboy problem will be: 



Cost 



CCS) 



C_(S-X), if 0 < X < S, and 

) 

c (X-S), if X > s, 



where X is the realization of demand for that period, and 
S is the inventory level in hand. Here S represents the 
decision variable for this problem. 
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The expected cost E[C(S)] which results from keeping 
an inventory level S in hand, will be: 



S 

E [C (S ) ] = Z C c (S-X)p n (X) + 2 C (X-S) p (X) (1) 

X=0 U X=S+l so u 



for discrete demand distributions and, 

S 

E [C (S) ] = C c / (S-x) f D (x) dx + C sQ / (x-S) f D (x)dx, (2) 
0 S 

for continuous demand distributions. 

Now, the expected cost for a period is to be minimized 
and we wish to find an optimal value of S, say S Q , that will 
yield the minimum expected cost. In other words, (1) and 
(2) are to be minimized with respect to S. 

In the discrete demand distribution case, two necessary 
conditions for a minimum at S Q are: 



E[C(S q + 1)] - E[C(S o )] > 0, and 
E[C(S q - 1)] - E[C(S q )] > 0. 

These conditions yield the well known result [7]: 
Pd^o' 1 ) - c c + C so lPD< S o>’ 



(3) 



as conditions for a minimum expected-cost solution at S Q . 

Similarly, for continuous demand distributions, to 
minimize (2), the first derivative of E[C(S)] is taken with 
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respect to S, and after it is set equal to zero, yields 



W 



c 

so 




so 



( 4 ) 



where Fp(S 0 ) is the cumulative distribution function (CDF) 
of demand at S Q . Thus, the S Q values which satisfy equations 
(3) and (4) are the optimum inventory levels to keep on hand 
in order to minimize the expected cost. 

This solution (under risk) is quite idealistic, since it 
assumes that the distribution of demand is known. In fact, 
in reality we are seldom certain about the distribution. In 
most cases, neither the distribution of demand, nor an ade- 
quate amount of data to provide a good basis for estimation 
of the distribution is initially available. The situation 
when: 

1. Distribution of demand is not known, and 

2. Enough data from the system to estimate the distri- 
bution of demand is not available, but 

3. We do have some idea about the range of demand D, 

say 0 < D < D , is one we shall call "uncertainty." 

max J 

B. SOLUTION OF NEWSBOY PROBLEM UNDER UNCERTAINTY 

It has been shown that whether minimax cost, minimax 
regret, or the Laplace approach, assuming a uniform distri- 
bution of demand is chosen as a principle of choice, the 
optimal decision under uncertainty is the same [6], [7], 

namely. 
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S o = 



C D 
so max 

C + C 
c so 



, for continuous demand. 



S < 
o 



C „(D + 1) 

so 1 max J 

C + C 
c so 



< S D H- 1, 



for a discrete demand distribution. The results of these 
uncertainty procedures are not affected by newly generated 
demand data unless it causes one to revise the range of demand. 

C. TRANSITION PHASE OF THE NEWSBOY PROBLEM 

Now the two bounds or limit solutions (uncertainty and 
risk conditions) for the transition phase of the Newsboy 
problem are reviewed. Implicitly, to proceed in accordance 
with these current solutions, one begins under Uncertainty 
and uses the Uncertainty optimum inventory level S Q (based 
on the range of demand) , until demand data is adequate or 
complete to estimate the distribution of demand and change 
to the optimum inventory level S Q for the Risk conditions. 

As mentioned before, the purpose of this study is to find a 
good transition phase procedure for the intermediate periods 
when some data has accumulated and is available. A transition 
phase procedure may be characterized by two decisions: 

1. When to begin the transition phase (i.e., at what 

period should we leave Uncertainty procedures to begin to 

use an estimator of S ) . 

o' 

2. What estimator of S to use. 

o 

After explaining the nature of the transition phase in 
the next chapter, three candidate estimators for the optimum 



12 



inventory level S Q to use during transition phase will be 
proposed. 

As we remember, the optimal choice of S under Risk was 
based on a quantile. Let us call, 



P = 



so 



C + C 
c so 



and then from equation (4) we have. 



S 



o 




f^Cp) , 



which provides a minimum expected-cost solution under Risk. 
Thus, the problem is to locate the p^ quantile of the 
Cumulative Distribution Function, F^X) . As more and more 
data is collected, one might be able to estimate, 



S 



o 



/\ 




F (p) , 



as a function of the data. Some alternative approaches to 
estimate this quantile were considered to be: 

1. To consider similar systems for which we have demand 
data, and then to apply a Bayesian approach. 

2. To decide about the form of the distribution of 
demand and carry out a maximum likelihood procedure to find 
the parameters . 

3. To apply a nonparametr ic approach. 

For the first two cases, if one is able to decide about 
the form of distribution such as lognormal, gamma, compound- 
poisson process, beta, etc., then the p t ^ 1 quantile of demand 
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might be estimated quite well subject to data. It is the 
third approach, that of nonparametr ic statistics, which is 
the center of interest in this thesis, and forms the basis 
for the chapter that follows. 
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III. A NONPARAMETRIC APPROACH FOR THE NEWSBOY 
PROBLEM, TRANSITION PHASE 

In this section the possibilities for applying a non- 
parametric or distribution- free approach to the problem of 
the transition phase of the Newsboy problem will be explored. 
After explaining the reasoning for a nonparametr ic approach, 
three candidate estimators of the optimum inventory level 
S Q (or X ) will be proposed. 

A. WHY A NONPARAMETRIC APPROACH? 

As mentioned at the end of the preceding section, if the 
decision maker is able to predict the form of the demand 
distribution, then a parametric approach can be more efficient 
and have a higher convergence rate to the true value of the 
estimator than a nonparametr ic equivalent. Hadley and 
Whitin [5], give various considerations relating to choice 
of a demand distribution and about inventory demand predic- 
tions, in detail. For a parametric approach, such as would 
result from fitting a theoretical distribution to data, they 
suggest to use poisson distribution when the lead time is 
relatively low, or normal distribution when, it is not rela- 
tively low, etc. 

Since it is quite difficult to specify the form of the 
demand distribution, a nonparametr ic approach was considered 
to be of particular interest for this problem. Also, Hadley 
and Whitin point out that, for single-period inventory models, 
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like the Newsboy problem, the empirical distribution function 
(nonparametr ic) is reasonably applicable, unless we are forced 
to operate the model at the tail of the distribution. [5]. 

Performance of the model decision rules was assumed to 
be critical during, say, the first 50 periods when data is 
inadequate. The decision maker needs to operate the system 
well during that critical time. Later on he will have a 
sufficient amount of data to treat the problem as one under 
Risk and simply use the optimal decision rules for Risk. For 
these early periods (transition phase) a nonparametr ic approach 
was used. 

B. THE EFFECT OF THE ESTIMATORS OF OPTIMUM 

INVENTORY LEVEL ON EXPECTED COST 

Before we suggest some estimators of optimal inventory 
level, X , it may be useful to see how those estimators 
affect the total expected- cost in a period. In the work 
that follows, it is assumed that the demand process is time 
independent (stationary) , which means that seasonal effects 
and trends are removed. Also, we assume that for the previous 
period we know what the demand was. 

Before we continue on with the mathematical development, 
it is useful to expand on the conditions for an optimal 
expected-cost solution. When the demand distribution is 
assumed known, we remember from equation (2) that the expected 
cost for the case of continuous demand is: 

S » 

E [C (S) ] ■ C c ! (S - x) f D (x)dx + C so / (x - S) f D (x)dx. 
o S 
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To find an optimum inventory level S q , taking 



E'[C(S)] 



d E |~C (S) ] 
d S 




( 5 ) 



and setting it equal to zero (necessary conditions for 
minimization) , we obtain: 



- C + C F n (S ) + C F n (S ) = 0 , or 

so so o' c D o' ’ 



W - 



SO 



C + C 
c so 



= p , and S Q = F^ 1 (p) . 



Similarly, 



E"[C (S) ] = d E l c ( s )] 



d S' 



S = S 



= C f n (S ) + C f n (S ) 
so o' c D o' 



= (C + C ) f n (SJ > 0 , 
v so c D o' 



( 6 ) 



which forms a sufficient condition for minimization. Here, 

4 - 1 - 

S Q is the value that corresponds to the p quantile of the 
population, say X . 

We are interested in the case where the population is 
not known, but there exists a sample of size n of demand 
data, acquired during the last n periods. 

Let X^, X n be realization of demand during the 

last n periods. Although we don't know the value of X^ yet, 
we do have a chance to estimate this value as a function of 
available data, say: 



x p = g(X 1 .X 2 



•>v 








( 7 ) 
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To express the expected cost as the function of X , let 

6 = S - S = X - X . 
o o p p 

Equation (2) then becomes: 



s + 6 

E[C(S q )] = C c /° (S Q + <5 - x) f D (x)dx + C S0 



* C so i ( Y S o' 3) f D (x)dx (8) 

b + o 
o 

If equation (8) is expanded using the Taylor series around 
S Q , it gives 

E[C(S q )] = E[C(S o )] + E " [C (S Q ) ] E[5] + ElC(S Q )] + 



Since E'[C(S q )] = 0 by equation (5), and having E~[C(S )] 
from equation (6) , then for three terms the Taylor expansion 
is : 

E[C(S 0 )] 1 E[C(S 0 )] ♦ CC so *C c ) f D (S o) £ifi (9) 

Equation (9) shows that: 

1. Under Uncertainty conditions the expected cost for 

t h 

a period is linear in the mean square error of the p 
quantile estimators, i.e., in 

E[6 2 ] = E[(X p - X p ) 2 ] 

2. If one is able to find a consistent estimator. 

A 

Xp(n), of X p , then when n gets large the mean square error 
will converge to zero and E[C(S q )] will converge to E[C(S )], 

A 

providing that the moments of X p (n) exist and behave properly. 
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3. The optimal performance for an estimator here is to 



yield the expected total cost under Risk, E[C(S q )]. 

Now, we wish to find a good consistent point estimator 

X for X . Also, it is not required, but it is desirable to 
P P 

obtain one that is an unbiased estimator of X^. If so, 

/s 2 a 

E[(Xp - X ) ] will be equal to the variance of estimator X^. 

Then, by equation (9), if it is possible to improve this 

estimator reducing its variance, it will directly reflect 

and will decrease the expected cost accrued by using that 

estimator . 

To obtain candidate nonparametr ic estimators of X , 
order statistics will be used. In the following section 
we shall develop three different extimators of X^ , using 
an order statistics approach. Additional details about order 
statistics and empirical distribution functions related to 
quantile estimation are given in Appendix A. 

1. The First Candidate Estimator of 



Optimum Inventory Level: X ^ r ^ 



Let X^ < X (2) < 



. < X, be a set of n reali- 
gn) 



zations from an unknown distribution in order statistics 

t h t h 

form, where X- ^ is the r largest value (r order 
statistic) . Let r be defined in the following way: 



r = 



np , when np is integer 
[n p + 1] ,. otherwise. 



where p is the quantile, n is the number of periods so far, 

and [X] denotes the integer part of X. For large n, X^ r -j 

is a consistant and unbiased point estimator of X^ (Appendix A). 
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In Figure 1, the r™ order statistic X^ is represented 
graphically with the distribution function F^(X) and the 
empirical distribution function G^(X) . 

Know results of this estimator include the following. 
For a uniform (0,1) distribution, 



E [X 




E [ U (r)l 



r 

n+1 



( 10 ) 



and 

c ov[X r _ v X r ] = cov[U r v U f -J = -• ( ~ n ~ S 2 1 ' ) , for r < s 

irj LSJ <- rj lSJ (n+l) Z (n+2) 

(ID 

Similarly, for any other distribution, 



E[X (r) ] 



O* IT 

'X, r 



D 






( 12 ) 



and 



covFX X 1 ^ n ~ 5 + 1 ) . 1 

frl , f s) ^ 2 “It _ i e 

t J t J (n+l) Z (n + 2) f„(F n 1 (-I T ))- f n (F„ h-fr)) 



D" D '-n+l- 
(13) 



D'- D '■n+l 

provided r < s. 

Later these results will be used to compare the 
efficiency of the estimators of X . 

2 . The Second Candidate Estimator 
of Optimum Inventory Level: 

As a second candidate estimator of X , X, is 

P (r) 

offered to smooth the empirical distribution function by a 
linear interpolation. It might be considered that the area 
between G^(X) and the actual distribution function F^(X) 
may be a measure of deviation. So, here interpolation was 
done at the midrange, yielding a new empirical distribution 
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function R, which would possibly have a smaller area between 

R and F^(X), than the area between Gp(X) and F^(X) . A possible 

realization for X^ r ^ is shown in Figure 2. Since the area 

between those distribution functions is subject to change 

because of randomness, the results for that case would be 

evaluated numerically rather than analytically. 

3 . The Third Candidate Estimator 
of Optimum Inventory Level: 

As a third candidate estimator of X , a pooled order 

P 

statistic, combining other order statistics, is proposed. 

For such an approach to be successful by equation (9) , it 

has to reduce the mean square error to a lower value than 

the r 1 "* 1 order statistic does, if it is possible to make this 

pooled estimator an unbiased estimator of X . Then the 

problem is one of comparing the variance of this estimator 

t h 

with the variance of the r order statistic. 

To examine this subject, it is reasonable to begin 
with the uniform (0,1) distribution, since the expectations, 
variances, etc. of its order statistics are simple and 
exact, and then to discuss the assumptions that are necessary 
before we generalize it for other distributions, 
a. Xj- r ^ for the Uniform Distribution 

Let X , X (r+1) the order statistics 

from the uniform (0,1) distribution. In addition to the r^* 1 
order statistic, a pooled estimator in the form of a linear 
combination of the (r-1), (r) , (r+1) order statistics will 
be considered. Here cur third candidate estimator is: 
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FIGURE 2. Second Candidate Estimator 



X (r) = X p " a X (r- 1 ) + 6 X (r) + Y X (r + 1 ) » 



where a + 3 + Y = 1 



As mentioned earlier, if this estimator can be an unbiased 
estimator of X , then a comparison of the variances of the 
r ^* 1 order statistic (First candidate estimator), and this 
pooled estimator (Third candidate) will indicate which one 
has better performance. 

To maintain unbiasedness, we must have: 



E[X (r) ] = E[X (r) ] , or 

E[f (r) ] = a E[X (r . 1D ] + 3 E[X (r) ] + Y E[X (r+1) ]. 



If demand has a uniform distribution, then 

E i f (r)i ■ “ * 6 Cnrr> + Y • 

*£• 

This quantity is equal to E[X^ r ^j = only when a = y • 

One way to decide on the values of a, 3 and y is to choose 
values so that the variance is minimized. So, we would like 
to minimize 

Var (f ( r)) = Var(a X^.^ + 3 X (r) + Y x (r+1) ) ’ 
subject to, 

a + 3 + Y = 1 and a = y 

By the constraints on a, 3 and y, these coefficients can be 
written in terms of a only, yielding 

Var(X (r) ) - Var{a(X (r _ 13 + X (r+1) ) + d-2a) X (r) } , 
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or, 

Var(f (r) ) = a^VarCX^.^) + Var(X (r+1) )) 

+ 2 Cov (X x ^ » x ( r+ i) ) 

+ (2a-4a 2 )(Cov(X (r _ ir X (r) ) + Cov (X (r) ,X (r+1) ) ) 
+ (1 - 4a + 4a 2 )Var(X ( - r . ) ) .. 

Taking the derivative of Var (X ) with respect to a and setting 
it equal to zero gives, 

a = 2 Var (X^ r -j ) - Cov(X (r _ 1 - ) >X (r )) ~ Cov ( X (- r -) ,X fr + 1 -)) 

Var (X ( r_i-j) + Var (X (r+1) ) + ZCovCX^.^ ,X (r+1) ) - 

4[Cov(X r _ I0 ,X (r) ) + Cov(X (r) ,X (r+1) ) - Var(X (r) )] 

(14) 

Having available the values of variance and covariances by 
equation (11), equation (14) becomes 

2r(n-r+l) - (r-l)(n-r+l) - r(n-r) 

a = 

(r-1) (n-r+2) + 4r(n-r+l) + (r+1) (n-r) - 4(r-l) (n-r+1) - 

4r£n-r) + 2(r-l)(n-r) 

the n,r terms in the numerator and in the denominator cancel, 
resulting in a = 1/2. So, to minimize the variance of X, . , 
the coefficients have to be a = y = 1/2, and 6=0. 

This says that when a population is uniform, a third 
reasonable estimator of Xp is: 

f (r) = ^(r-D + X (r*l/ ’ (15) 
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yielding 



VaT < x (r) 



2r(n-r+l) - (n+1 ) 

2 (n+1) 2 (n+2) 



If the ratio of the two candidate estimators is taken. 



Var(X f O 2r (n-r+1) 

_iU — = > 1 , 

Var (X £ r ^ ) 2r (n-r+1 ) - (n+1) 

always. This implies that for a uniform distribution, 

Var (X^- r -j) > Var (X^), and is the better estimator of 

Xp. In Figure 3, for n - 20, these variances are sketched. 
When r is close to bounds (either close to 1 or n) , the 
difference between those two variances become larger. 




FIGURE 3. Variances of X ( . and X, . for Unif orm( 0 , 1) 
Distribution ,wnen n is Twenty . 
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These results apply, of course, only to uniform 
demand distribution. In the following paragraphs the gener- 
alizations to any demand distribution will be discussed. 

b. Generalization of to any Demand Distribution 

For distributions other than uniform, it is 
difficult to claim that X^ r ^ , as given by equation (15), is 
exactly unbiased. Also, the variance and covariances, as 
equation (13), will be functions of the density f^(X). Even 
though these two problems exist, some reasonable approxima- 
tions can be discussed. If we have both a large n and fairly 
flat densities, since (r-1), (r) , (r+1) are the closest 

neighbors among the other order statistic values, the densi- 
ties corresponding to these points would tend to equal each 
other and during the calculation of a, as in equation (14) , 
would essentially cancel. 

If we have flat densities (where the three points 
are nearly linear) , or if we have a large number of observa- 
tions (when the points are close together accordingly), 

Xj- r -j would be approximately unbiased, too. Thus, it may be 
reasonable in some cases to use X , . for other distributions 
besides the uniform. 

Incidentally, we would not like to try another 
pooled estimator of X^ such as 



(r) 



= a. 



X (r=2) +a l 



X (r-l) + f3 



X, , + 
(r) 



Y 1 A (r + 1) + Y 2 X (r+2) 



For this estimator, the assumptions which were made to permit 
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X(r) to 'apply the distributions other than the uniform may 
no .longer be appropriate. For example, in equation (13), 
f D (Fq 1 (j^y) ) and f^F^ 1 (£py) ) would probably have quite 
, different values, because these are the densities of the two 
farthest order statistics, X( r _ 2 ) aa d ^( r +2) 

Otherwise, if these densities were known, we would not use 
a nonparametr ic approach. 

So far, solutions of the Newsboy problem under 
"Uncertainty" and under "Risk" have been reviewed, and 
possible procedures for a transition phase between them 
discussed. Also in this chapter, three candidate estimators 
of the optimal inventory level S Q (or X^) were developed for 
use as alternative transition phase procedures. In the next 
chapter, the performance of these alternative transition 
procedures will be evaluated by a computer simulation of 
Newsboy problem decision making. Finally, from the simula- 
tion we will obtain information about what time to leave the 
minimax (Uncertainty) rule, for each transition phase pro- 
cedure . 
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IV. EVALUATION OF THE CANDIDATE ESTIMATORS 
OF OPTIMUM INVENTORY LEVEL, FOR THE 
TRANSITION PHASE OF THE NEWSBOY PROBLEM 



In Chapter III, the rationale for using a nonparametr ic 
approach for the transition phase of the Newsboy problem was 
discussed. It was shown that under Uncertainty the expected 
cost for a period was approximately linear in the mean square 
error of the optimum inventory level S Q (or X ) corresponding 

■f* V 

to the p quantile, and three candidate estimators of or 

S were introduced, 
o 

In this chapter, the performance of these estimators 
will be tested for several demand distributions by simulation. 
After explaining our reasons for using simulation, we will 
explain the structure of the simulation program used. In 
part C, the simulation parameters used in the program are 
described, together with the reasons for choosing these 
parameters. Part D illustrates the results of simulation. 
These results will be used later to compare the performance 
of the estimators of X^ and the transition phase procedures, 
and to see whether there were any apparent adverse effects 
of approximations or assumptions. 



A. RATIONALE FOR SIMULATION 

The purpose of the study was to find a transition phase 
procedure (i.e., which estimator of X^ or S Q to use and when) 
during the transition from Uncertainty to conditions under 
Risk for the Newsboy problem. As explained before, the 
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performance of an alternative transition phase procedure 
depends directly on the performance of the candidate esti- 
mator of X used, in that there is a one-to-one correspondence. 
P 

A primary reason for simulation, of course, was to study 
the question of how long one should continue using the mini- 
max procedure (for Uncertainty) before changing to a transi- 
tion phase using a quantile estimate. Simulation over several 
demand distributions should provide insight regarding when 
(with respect to cost) this change in the Newsboy problem 
should occur. Additionally, we wish to further examine the 
question of which of the candidate estimators to use. 

However, it was hoped that .the estimators of or S Q 
would come close to having desired order statistics proper- 
ties, being exact for the case of the uniform distribution 
(equations (9), (10), (11) and (12)). For example, for the 

r £ ^ order statistic X, . it has been shown that 

(r) 

F ri (E[X, .1) = — 7 -T for a uniform distribution, but 
D (r) n+ 1 

F D^ E f X (r)^ - iiTT £or convex > and F D^ E f X (r)^ - nTl for 
concave cumulative distribution functions of demand. [4]. 

Unfortunately, there is not an analytical measure to show 

how much the slack in these inequalities would change 

depending upon the amount of convexity or concavity of the 

demand distributions. 

Also, for the third candidate estimator of X , some 

P 

additional assumptions were made that might not hold for 
every distribution. 
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In summary, our reasons for simulation are: 

1. To find the proper time (period) to begin the 
transition phase (to use one of the candidate estimators) , 
leaving the "Uncertainty" case and the minimax decision rule, 

2. To provide information to supplement incomplete 
analytical knowledge of order statistics regarding popula- 
tion quantiles, 

3. To evaluate the impact of the approximations made 
in developing the third candidate estimator, X^-^, over 
various demand distributions. 

B. DESIGN FOR SIMULATION 

The simulation program was written in the FORTRAN IV 

language. In every simulation run, inventory stocking 

decisions are made successively using one particular demand 

distribution and one particular value of the quantile 
C so 

(P = — p — — ) , say, demands generated from a normal parent 
L c so 

for quantile 0.3 ... etc. The demand distributions and 
various quantile values used in the simulation are given in 
the next section. Each run consisted of 50 successive 
decision periods starting with no demand information other 
than the range estimate, and replicated 60 times. 

As output values of the program, costs for each period 
were averaged over the number of replications. These average 
costs per period were calculated for five different cases: 

1. The Uncertainty case, with the minimax decision rule 
being used, denoted by M. 

2. The ideal case (Risk) , with the quantile value known, 
denoted by ID . 
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3. The transition phase case using estimator X 
by R . 




4. The_transition phase case using estimator X 
by R. 



( , denoted 

(rj 



The minimax (for Uncertainty) and ideal (for Risk) cases 
would give constant inventory levels, i.e., the same amounts 
all the way through the simulation run. For Uncertainty, the 
range of demand was approximately known and would not change 
unless it was revised, while for Risk conditions it was 
assumed that the distribution of demand, and thus the quan- 
tile, was known. 

Transition phase optimum inventory level values generated 
by ^C r ) ’ X (r) an< ^ X (r) were subject to change from period to 
period until finally they converge to the value associated 
with the Ideal case. For the first two periods, X^-j, X^.^, 
X(r) would not lead to an inventory level, since they need 
at least two data points. 

The program was run for SO periods, and replicated 60 
times. So, period costs for each transition phase procedure 
were averaged over 60 replications. Program steps can be 
ordered in the following way: 

1. Define the demand distribution yielding the quantile 
p, give the stockout cost C SQ , the carrying cost C c , and the 
range of demand. 

2. Compute and set minimax and ideal case optimum 
inventory values. 

3. Generate a random demand value from the specified 
demand distribution. 
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4. Find the single-period cost for each case (minimax, 

, X (r) and Ideal) by comparing their inventory 

levels and the demand generated for this period. Record 

these costs. Find the next period inventory levels revising 

X/- x"/- X"/- and go to the next period. 

(r) » (r) (r) & K 

5. Repeat steps 3 and 4 for 50 periods. 

6. Again, beginning from the first period, repeat 
step 5 for 60 times (replications) , and at the end of 60 
replications find the average costs, and the sample variances 
of costs over 60 replications, for each period. 

7. Tabulate the results. 

The results obtained from this simulation will be examined 
in part D. 



C. SIMULATION PARAMETERS 

In this simulation a variety of possible forms of distri- 
butions, including continuous, discrete, symmetric, skewed 
(long-tailed), etc., were chosen. Demand distributions and 
their parameters used in this simulation were the following: 



Distribution 

Uniform 

Poisson 

Chi-square 

Lognormal 

Normal 



Parameters 

the range over [0,30] 

with mean 10 units 

with 10 degrees of freedom 

generated by a normal parent 
having a mean of 3 and a variance 
of 1 

with mean 35 and variance 100 
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For D in the Minimax procedure, \i + 3 x is used for 
max _ 

'-'SO 

tailed distributions. Quantiles p = ^ ^ — , tested were 

c So 

0.10, 0.30, 0.50, 0.70, 0.90, 0.95 for each distribution. 

These were determined so that C SQ + C c was constant (in this 

simulation C + C = 4.0 was used). So, for each (quantile, 
so C ~i j 

distribution} combination the differences of the expected 
period costs would come from the performance of the estimators 
of Xp, as was shown earlier by equation (9) . 

In the next section we shall present and discuss some 
of the results from the simulation. 

D. SIMULATION RESULTS 

In this part, simulation results will be evaluated and, 
for the Newsboy problem, a transition phase procedure will 
be proposed. In presenting these results we will identify 
the quantile value as p. Tables I and II give overall mean 
costs and standard deviations of cost. Then, Table III will 
show at what period the transition procedure should begin 
for a particular distribution, and quantile. Table IV will 
give relative efficiency of the best transition procedure 
to the Ideal case (decision under Risk) . 

On the next page. Table I gives the overall mean values 
of the cost per period (over 50 periods and 60 replications) 
for each transition case, using C sQ + C c = 4.00. As can be 
seen from Table I, there does not seem a significant differ- 
ence using R, R, or R, even though R looks quite consistent 
at the tail ( p= 0.95). One reason may be the fact that the 
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TABLE I 



Average Costs Over 50 Periods and 60 Replications 

having C + C =4.0 
& c so 

M=Minimax, R=X^ r -j, R=Xj- r -j, R=X^ r ^, ID=Ideal Cases 



p 




Uniform 


Poisson 


Chi-Square 


Lognormal 


Normal 




M 


5.49 


3 . 23 


3.01 


67.00 


11.66 




R 


6.35 


2.38 


2.91 


14.50 


8.09 


0.1 


R 


6.31 


2.38 


2.88 


14.50 


8 . 09 




R 


6.30 


2.38 


2.87 


14.48 


8.07 




ID 


5.49 


2.04 


2 .56 


11.44 


6.93 




M 


12. 68 


4.60 


5.92 


255.96 


20.04 




R 


13.58 


4.55 


5.98 


36.34 


14 .25 


0.3 


R 


13.55 


4.52 


5.96 


36.23 


14.18 




E 


13.46 


4.56 


5.94 


36.38 


14.16 




ID 


12.68 


4.26 


5.62 


30.15 


13.17 




M 


15.08 


5.97 


11 . 08 


339.86 


15.71 




R 


16. 03 


5.35 


7.47 


51.22 


15.80 


0.5 


R 


15.97 


5.33; 


7.45 


51.37 


15.78 




R 


15 . 91 


5.34 


7.45 


51.35 


15.73 




ID 


15.08 


4.98 


6.98 


43.49 


14.80 




M 


12.75 


7.33 


11.81 


298.67 


14.97 




R 


13.75 


4.74 


7.21 


57.24 


14.24 


0.7 


R 


13.66 


4.74 


7.21 


57.23 


14.19 




E 


13.60 


4.72 


7.20 


57.66 


14.19 




ID 


12.75 


4.44 


6.67 


48.72 


13.24 




M 


5.39 


4.37 


6.26 


132.42 


9.33 




R 


6.31 


2.68 


4.40 


42.36 


7.76 


0.9 


R 


6.30 


2.67 


4.38 


42.38 


7.76 




R 


6.30 


2.66 


4.36 


42.58 


7.73 




ID 


5.39 


2.33 


3.71 


35.91 


6.70 




M 


2.83 


2.38 


3.34 


70.42 


5.27 




R 


3.77 


1.75 


3.00 


31.80 


5.04 


0.95 


R 


3.76 


1.75 


2.99 


31.74 


5.02 




E 


3.76 


1.75 


3.00 


32.44 


5.04 




ID 


2.83 


1.41 


2.27 


25.47 


3.91 
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difference between variance of and variance of X^ , 

was not much (Figure 3). But if one wishes to see the dif- 
ference no matter how much it is, the following normalization 
by ideal costs may be considered: 

1. For a distribution from Table I, for every quantile 
R.'-ID 

£ find — j-g — , and sum this value over all quantiles for 
every R^, where R^ is R, R or R. 

2. Repeat step 1 for every distribution. 

3. Pick the R^ ' s which have the smallest sum among 
R, R, R for every distribution as their best procedure. 

The results of this normalization were as follows: 





Uniform 


Poisson 


Chi-square 


Lognormal 


Normal 


R 


0.871 


0.767 


0.859 


1.253 


0.839 


R 


0.845 


0.752 


0.831 


1. 251 


0.824 


I 


0.828 


0.755 


0.821 


1.295 


0.816 



The best transition procedures were to use R for uniform, 
R for poisson, R for chi-square, R for lognormal, and R for 
normal distributions. To see the effects of the very right 
tail of a distribution on transition procedures, p = 0.99 
was used for normal and chi-square distributions. At this 
quantile the minimax costs were below R, R, and R costs. 

This says minimax can be used at the very tail of a distri- 
bution along the transition phase, providing very accurate 
prediction of the range of demand (say +5%). Otherwise, if 
we remember from Table I for p = 0.95, R, R and R were 
better than minimax approach. 
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From Table II (overall standard deviations over 50 periods 
and 60 replications) , it is seen again that there is not much 
difference of standard deviations, even though R has a better 
trend (smaller sample variance of cost) among .R, R and R. 

From these two tables it can be said that there is not 
much difference of the costs yielded using R, R and R, but 
if an order is necessary it should be the one using R for p 
less than 0.95, and R at the tail (i.e., when p is greater 
than 0.95) . 

During the simulation, as the number of periods increase 
the estimators converge to X^ or the true S Q , while the 
minimax and ideal case inventory levels were constant. 

Since the values of X^, X^ , X^ usually change from 
period to period because they are updated by the most recent 
demand value, they will generally be associated with a larger 
variance of cost than the minimix and ideal cases. This is 
reflected in Table II. Sometimes the minimax procedure had 
a smaller sample variance of cost (while possibly yielding 
the largest average cost) than the ideal case. This may be 
explained as follows: 

The ideal solution we are using relates to minimizing 
expected cost, not the variance of cost. These two minimi- 
zation points are generally different, depending on the 
distribution function. If the minimax solution is different 
from the ideal (expected cost) inventory level and closer to 
the minimum variance solution, it would yield a smaller 
variance of cost. Further derivations about the variance 
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TABLE II 



Sample Standard Deviations over 50 Periods 
and 60 Replications 

M=Minimax, R=X^ r -j, R=X^ r -), R=Xj- r -j, ID=Ideal Cases 



p 




Uniform 


Poisson 


Chi-Square 


Lognormal 


Normal 




M 


3.15 


1.26 


1.79 


43.46 


3.77 




R 


6.36 


2.69 


3.27 


27.46 


9.89 


0.1 


R 


6. 26 


2.65 


3.21 


27.47 


9.70 




R 


6.26 


2.59 


3.18 


27.47 


9.53 




ID 


3.15 


1.62 


2.26 


15.90 


6.62 




M 


7.46 


3.18 


4.61 


76.10 


10.33 




R 


9. 00 


3.60 


4.93 


57.64 


12.17 


0.3 


R 


8.84 


3.53 


4.89 


57.69 


11.97 




R 


8.85 


3.51 


4.86 


57.71 


11.99 




ID 


7.46 


3.16 


4.42 


45.93 


10.95 




M 


8.82 


4.30 


6.35 


66.95 


12.34 




R 


10.31 


4.14 


6.37 


84.29 


12.73 


0.5 


R 


10. 23 


4.12 


6.33 


83 . 91 


12. 68 




R 


10.23 


4.14 


6.44 


84.23 


12.75 




ID 


8.82 


3.84 


5.80 


72.72 


12.09 




M 


7.30 


3.47 


4.87 


46.79 


9. 62 




R 


9.08 


3.95 


6.66 


99.68 


11.65 


0.7 


R 


9.01 


4.04 


6.83 


99.68 


11.65 




R 


8.87 


4.09 


6.76 


99.23 


11.67 




ID 


' 7.30 


3.55 


5.66 


92.87 


10.42 




M 


3.16 


1.25 


1.87 


32.45 


3.90 




R 


6.59 


3.28 


5.72 


99.99 


9.10 


0.9 


R 


6. 55 


3.33 


5.76 


99.23 


9.09 




R 


6. 55 


3.29 


5.69 


98.05 


9.01 




ID 


3.19 


2.57 


4.10 


94.42 


6.34 




M 


1.66 


0.69 


1.12 


31.66 


2.07 




R 


6.50 


2.91 


5.36 


94.64 


8.50 


0.95 


R 


6. 50 


2.95 


5.34 


92.85 


8.47 




R 


6. 50 


2.94 


5.23 


92.74 


8.46 




ID 


1.66 


1.49 


3.15 


84.97 


4.11 
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minimization solutions of the Newsboy problem are given in 
Appendix B. 

The simulation results about the period to begin a 
transition phase are given in Table III. Shown are the 
first periods in which the average period cost (over 60 
replications) using the estimator was below that obtained 
using the minimax rule. In the simulation it was observed 
that once costs using an estimator fell below the minimax 
costs, they stayed below for the remaining time periods. 
Simulation results showed that, on the average, at the fifth 
period we were supposed to begin a transition procedure, 
ceasing to use the minimax procedure. In other words, for 
the first five periods the large variability of candidate 




the minimax procedure yielded. As a very special case, if 
the distribution of demand is uniform, minimax and Ideal 
case solutions would match each other and thus the uniform 
distribution is omitted from Table III. Also by coincidence, 

^SQ 

if p = p — 7Tr — is at ° r close to an intersection of the 
L c L so 

implied uniform CDF of minimax procedure and the actual 
distribution function of demand, then the minimax procedure 
may yield the Ideal case solution. Otherwise, minimax in- 
ventory level has always a higher cost than the Ideal case. 

In Table III the large numbers are those specific values 
which relate to a quantile P that is close to minimax- - 
actual CDF intersection point. So it takes time for R, R 
and R to yield costs less than the minimax. No matter how 
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TABLE III 





Periods to Begin Transition Phase 
Leaving the Minimax Procedure 


Quantile 


(p) Poisson Chi-Square Lognormal Normal 


0.1 


6 6.5 6 


0.3 


10* 10* 4 4 


0.5 


6 4 4 10* 


0.7 


4 4 4 6 


0.9 


5 5 4 5 


0.95 


6 5 4 6 



Ideal -Minimax intercept points (or, close to these points) 



40 



long it will take to fall below the minimax cost (here, 
nearly at period ten) from simulation outputs it was observed 
that after the fifth or sixth period the cost difference 
between the minimax and the transition procedures was not 
large. This also can be predicted from the pattern values 
of Table III. The numbers in Table III are relevant for all 
candidate estimators, i.e., R, R and R all fell below the 
minimax cost simultaneously at the same period. 

As a measure of efficiency for a transition procedure, 
the ratio of average cost induced by this procedure to Ideal 
average cost is used. Table IV gives these ratios for the 
average costs over second and fifth 10 periods for R (since 
R seems the promising transition procedure among the other 
ones) . 

We will examine the' data of Table IV for the second 10 
period, assuming that the minimax has been left and the 
transition phase has begun. We also believe that convergence 
toward the true value of Ideal inventory level S Q for the 
third, fourth, ... 10 periods will not be worse than the 
convergence accomplished during the second 10 period. 

If the trans ition .procedure (here, R) was operated at 
the tail (i.e., p 0.9), then the efficiency ratio for the 
second 10 periods was on the average 1.16 and more (worse). 

As the quantile gets close to 1.0, the minimax procedure was 
yielding much better results than a transition procedure 
(assuming the range of demand is predicted quite well) . 
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TABLE IV 



The Ratio of the Second and Fifth Ten-Period 
(11 to 20, 41 to 40) Average Costs £o the 
Ideal Case (Under Risk) Costs, Using X, 

Transition Procedure ' 



Quantile 


Cp) 


Uniform 


Poisson 


Chi-Square 


Lognormal 


Normal 


0.1 


s 

F 


1.063 

1.000 


1.065 

1.023 


1.058 

1.046 


1.033 
1. 002 


1.108 

1.008 




S 


1.065 


1.059 


1.050 


1.040 


1.066 


0.3 




F 


1.007 


1.041 


1.003 


1.004 


1.008 




S 


1.056 


1.069 


1.039 


1.041 


1.077 


0.5 




F 


1.013 


1.021 


1.013 


1.004 


1.007 




S 


1.073 


1.069 


1.061 


1.057 


1.054 


0.7 




F 


1.000 


1.001 


1.006 


1.006 


1.004 


0.9 


S 


1.102 


1.113 


1.128 


1.190 


1.096 




F 


1.008 


1.016 


1.052 


1.026 


1.027 




S 


1.170 


1.193 


1.211 


1.240 


1.183 


0.95 




F 


1. 002 


1.040 


1.054 


1.062 


1.055 



Notes: . 

S stands for "SECOND, M F stands for "FIFTH." 

When p is less than 0.9 overall ratio is 1.060 for second 
and 1.010 for fifth 10 periods. 
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For the rest of the quantiles where p is less than 0.9, 

this ratio for the transition procedure was almost constant 

1.060. To test the effects of quantiles and the effects of 

the form of demand distributions on the efficiency ratio 

(1.060), a two-way Analysis of Variance could be applied to 

the data of Table IV for quantiles less than 0.9. 

Different distributions were significant at a level 0.03, 

and the different quantiles were significant at a level of 

0.01. This result was explained that for C + C =4.0 

r so c 

cost units, whatever the distribution of demand was and for 
any quantile less than 0.9, the ratio of the cost induced by 
transition phase procedure (X^-j , here) to the Ideal cost 
(expected cost under Risk) would be practically constant at 
1.060 for second 10 period, and would not exceed this amount 
for the further periods (for example, at the fifth 10 period 
this ratio was 1.01). 

On the next part, as the summary of this chapter, and 
using the results obtained in it, a Single-Period inventory 
(Newsboy) problem transition phase procedure- is suggested. 

E. PROPOSED NEWSBOY TRANSITION PROCEDURE ALGORITHM 

/ 

(1) Approximate the bounds of Demand (Range). 

C s o 

(2) Set p = f f — . 

c c + c s0 

(3) If p > 0.90, go to step 6, otherwise go to step 4. 

(4) Use minimax for the first 5 or 6 periods. 

(5) For 6 _< n _< 50, use R (unless a reason for R, R) , 
and go to step 10. 

(6) If p £ 0.95, go to step 8, otherwise go to step 7. 
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(7) If step 1 is accurage, may use minimax for all periods 
and go to step 10, but if it is not accurate, go to 
step 8. 

(8) Use minimax for the first 5 or 6 periods, then go to 
step 9. 

(9) For 6 <_ n £ 50, use R and go to 10. 

(10) Analyze data. 

*May fit a probability distribution to the data and begin 
Decision Under Risk case. 



*May use R, R or R for n > 50 too. Since they converge 
X quite fast, the risk of using them after 50 periods hope- 
fully will be nearly equal to the risk of Ideal case. 
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V. CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER STUDY 



For the Newsboy problem, a nonparametric decision pro- 
cedure that allows transition from Uncertainty to decision 
under Risk has been developed. Three estimators of the 

demand corresponding to the p™ quantile were obtained, 
c so 

where p = p — — p — . Testing these estimators over several 
L c + L so 

distributions of demand with Newsboy cost criteria, it was 
found that excluding cases such as drastic skewness, pooling 
the (r-l)st and (r+l)st order statistics for R provided 
essentially the best estimator, where R and R also worked 
well . 

It was also found that for the distributions tested and 

C 

excluding ( r - s ° r ) values very close to 0.0 or 1.0, that 
L so L c 

the switching Newsboy decision making from minimax solutions 
to quantile estimator could be made, possibly at the sixth 
period. 

For further work, some parametric studies such as log- 
normal distribution, or beta and gamma families, depending 
upon whether the population of demand is finite or not, may 
be tried. Then, given that the decision maker does not know 
what the distribution of demand is, which procedure, nonpara- 
metric or parametric, gives better response to Uncertainty 
may be compared. As another recommendation, some work can 
be afforded to improve nonparametric approaches further. 
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